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ABSTRACT 


Deep  water  directional  wave  spectra,  measured  by  an  NDBC  3-meter  buoy  off 
Montery  bay,  are  transformed  to  shallow  water  using  a  linear  refraction  model.  The 
transformed  directional  spectra  are  compared  with  measured  spectra  using  pressure 
gauge  arrays  in  shallow  waters  at  Marina  and  Santa  Cruz. 

The  classical  Longuet-Higgins  et  al.  (1963)  method  of  computing  directional  wave 
spectra  and  a  new  exact  Fourier  coefficients  representation  method  (Grauzinis,  1989)  are 
used  to  compute  directional  wave  spectra.  The  new  method  of  computing  directional 
wave  spectra,  which  represents  bi-modal  distributions  of  wave  energy  exactly  matching 
the  measured  Fourier  coefficients  to  second  order,  demonstrates  improved  directional 
resolution  over  the  classical  technique. 

This  work  examines  the  accuracy  and  limitations  of  modeling  linear  refraction  by 
comparing  with  field  observations  over  complex  bathymetry.  In  general,  linear  refrac¬ 
tion  can  give  reasonable  energy  and  direction  estimates  starting  with  deep  water  spectra, 
but  notable  exceptions  can  occur.  The  largest  prediction  error  occured  at  Marina  on 
18  January  1988  for  the  case  of  a  severe  storm.  This  is  presumably  due  to  diffractive  and 
non-linear  effects  of  the  high  waves  causing  loss  of  accuracy.  The  linear  refraction 
model  is  not  suitable  for  handling  such  problems. 
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I.  INTRODUCTION 


A.  MOTIVATION  OF  STUDY 

Wave  climatology  plays  an  important  role  in  coastal  dynamics.  A  basic  problem  is 
how  to  use  deep  water  wave  data  to  predict  shallow  water  wave  characteristics,  which 
would  allow  predictions  of  waves  at  any  nearshore  location  using  only  the  measurement 
at  one  deep  ocean  location.  But  there  are  two  associated  questions.  First  ,  how  to  char¬ 
acterize  ocean  wave  behavior?  Second,  how  to  transform  waves  from  deep  water  to 
shallow  water?  Since  directional  characteristics  of  ocean  waves  are  important  to  under¬ 
stand  for  wave  generation  and  wave  refraction  behavior,  we  choose  to  use  the  directional 
wave  spectral  density  (DSD)  for  characterizing  the  ocean  wave  behavior  and  we  will  use 
a  linear  refraction  model  to  transform  the  deep  water  DSD  to  shallow  water. 

The  area  of  Monterey  Bay  has  been  selected  for  this  comparative  study  of  ocean 
waves.  This  area  presents  an  interesting  problem  for  predicting  directional  wave  spectra 
because  of  its  variable  shoreline  orientation  relative  to  prevalent  eastward  waves  and  the 
presence  of  the  Monterey  submarine  canyon  which  is  the  largest  in  the  western  hemi¬ 
sphere  and  causes  substantial  perturbations  on  the  nearshore  wave  field.  Figure  1  shows 
the  irregrlar  bathymetry  of  the  Monterey  submarine  canyon  and  three  sensors  locations 
used  in  the  comparison  study. 

A  NOAA,  3m  buoy  is  used  in  deep  water  to  measure  the  three  orthogonal  compo¬ 
nents  of  the  wave  field  heave,  pitch,  and  roll  directly.  In  shallow  water,  square  slope 
arrays,  6m  on  a  side  composed  of  4  pressure  sensors  each  are  used  to  measure  and 
compute  the  equivalent  three  orthogonal  wave  components  at  Marina  and  Santa  Cruz. 

B,  HISTORICAL  BACKGROUND 

Historically,  Monterey  Bay  was  the  first,  and  presumably  one  of  the  most  intensively 
studied,  locations  for  wave  refraction.  Wave  refraction  diagrams  of  Monterey  Bay  were 
first  constructed  in  1948  by  Johnson,  O  Brien  and  Isaac  (Wiegel,  1964).  The  first  in¬ 
tensive  refraction  studies  were  performed  for  amphibious  landing  exercises  at  Fort  Ord 
(1954)  followed  by  studies  by  the  Corps  of  Engineers  (1958)  for  the  small  craft  harbor 
at  Santa  Cruz.  Wave  refraction  and  littoral  drift  calculations  have  been  performed  by- 
Dorman  (1969)  and  Amal  et  al.  (1973).  All  these  studies  were  performed  using  hand- 
drawn  refraction  diagrams.  More  recently,  numerical  refraction  calculations  have  been 
made  for  a  number  of  specific  marine  construction  sites  (e.g.,  Thornton  (1978,  1980)) 
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and  erosion  studies  (Oradiwe,  1986;  McGee,  1986).  The  present  work  is  an  extention 
of  the  work  by  Khalid  (1989)  who  used  the  low  resolution  spectral  estimator  method  of 
Longuet-Higgins,  Cartwright  and  Smith  (1963). 

C.  OBJECTIVE 

The  objective  of  this  study  is  to  investigate  whether  linear  refraction  can  be  used  to 
transform  directional  wave  spectra  from  deep  to  shallow  water  over  complex 
bathymetry.  Several  different  ways  of  computing  directional  spectral  density  (DSD)  are 
presented.  The  deep  water  DSD  are  transformed  to  shallow  water  using  the  refraction 
approach  originally  proposed  by  Longuet-Higgins  (1957).  These  transformed  spectra  are 
compared  with  DSD  measured  directly  in  shallow  water  for  which  the  data  had  been 
acquired  from  pressure  gauge  arrays  located  at  Marina  and  Santa  Cruz.  In  Chapter  II. 
of  this  thesis,  linear  wave  theory,  two  methods  of  computing  directional  wave  spectrum, 
and  wave  spectrum  transformation  using  the  theoretical  approach  will  be  addressed. 
Chapter  III.  describes  the  data  acquisition  and  sensors.  Chapter  IV.  contains  the  results 
and  discussion.  Conclusions  are  presented  in  Chapter  V. 


3 


II.  THEORETICAL  BACKGROUND 

A.  THE  LINEAR  WAVE  THEORY 
1.  Formulation 

Linear  surface  gravity  wave  theory’  has  been  successful  in  describing  many  ob¬ 
served  wave  phenomena.  In  the  formulation  of  this  theory,  see  for  instance  Kinsman 
(1965),  the  basic  assumptions  can  be  summarized  as 

•  wave  frequencies,  co,  are  much  greater  than  the  earth's  rotation  frequency,  Q 

•  the  fluid  is  homogeneous,  or  co>V,  the  Brunt-  Vaisala  frequency 

•  surface  tension  effects  are  negligible,  or  uc,Pa)  ,  where  is  the  frequency  of  the 
capillary’  waves 

•  the  fluid  motion  is  irrotational,  or  V  x  V  =  0,  which  allows  the  velocity  potential, 
<t>,  to  be  defined  in  terms  of  the  water  particle  velocity  vector  as  V  =  —  V(f> 

•  the  fluid  is  incompressible,  or  V  •  V  =  0 

The  combination  of  the  assumptions  of  irrotationality  and  incompressibility 
leads  to  the  basic  equation  of  linear  wave  theory,  Laplace's  equation 

V20  =0  (1) 


Analytical  solutions  to  this  equation  are  obtained  by  applying  the  following  boundary 
conditions: 

•  waves  are  periodic  in  space  and  time 

•  the  bottom  boundary  is  horizontal  and  impermable,  which  is  prescribed  by 

»U— : jj-U-o  P) 


•  the  linearized  kinematic  free  surface  boundary  condition  is  expressed  as 

4=- - 1,-0 

ct 


(3) 


and  the  linearized  dynamic  boundary  condition  as 

d<t>  , 

where  w  is  the  scalar  vertical  velocity,  h  is  the  local  water  depth  ,  tj  is  the  instanta¬ 
neous  free  surface  elevation,  g  is  the  acceleration  of  gravity,  and  i  is  time. 
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In  the  linearization  process,  it  is  assumed  that  the  wave  amplitude  a  is  small  when 
compared  to  the  wavelength  L  and  depth,  or  that  ~  <^1  and  -^-<^1.  These  restrictions 
make  the  linear  wave  theory  strictly  applicable  to  small  amplitude  waves  of  infinitesimal 
steepness.  The  resulting  solution  to  the  wave  equation  is  a  single  harmonic  function. 
Invoking  the  superposition  principle,  the  instantaneous  sea  surface  elevation  can  be  re¬ 
presented  as  an  infinite  sum  of  sinusoids  as  expressed  by 


QO  OQ 

n(*)  =  X'7" =  Tj°n  cos(Kn  •x-(o„t  +  en)  (5) 

n=l  «=l 

w'here  Kn  is  the  horizontal  vector  wavenumber,  x  =  (jc,^)  is  the  position  vector,  w„  is  the 
circular  frequency  and  cn  is  the  initial  phase  angle.  The  dispersion  relationship  is  obtained 
from  the  solution  of  the  wave  equation  in  applying  the  surface  boundary  conditions 

oj 2  =  gk  tanh(M)  (6) 

The  spectral  component  index  n  is  now  dropped  for  convenience. 

Two  limiting  regions  can  be  considered,  deep  and  shallow  water  wave  approxi¬ 
mations.  These  simplified  solutions  correspond  to  the  consideration  of  asymptotic  forms 

for  the  hyperbolic  functions  entering  the  basic  solution  definitions.  These  extreme  cases 

h 

are  generally  identified  by  the  ratio  —  ,  w’here  h  is  the  local  water  depth  and  L  is  the 
wavelengh  as  computed  by  the  expressions  corresponding  to  each  one  of  the  cases.  For 
h/L  >  1/2  the  solution  is  considered  as  belonging  to  the  deep  water  limit  and  for  h/L 
<  1/20  the  solution  corresponds  to  the  shallow  water  case.  In  between  the  full  solution 
must  be  considered  and  the  solution  is  known  as  the  intermediate  water  wave  solution. 

The  solution  of  the  linear  wave  equation  can  be  used  to  describe  the  wave  in¬ 
duced  motions.  In  particular,  the  relation  between  wave  elevation  rj  and  subsurface 
pressure  head  p  is  described  by 

]w+;]  (7) 

Open  ocean  data  are  invariably  collected  in  the  form  of  time  series  of  processes 
that  are  inherently  random.  Spectral  analysis  techniques  are  a  natural  first  choice  to 
study  these  data.  The  previously  presented  solution  to  the  Laplace's  equation  gives  rise 
to  simple  spectral  relationships  between  the  sea  surface  elevation  spectrum  and  pressure 
spectra  when  considering  the  water  column  as  a  constant  parameter  linear  system.  Un- 
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der  such  a  system,  any  two  qua.  cies  are  connected  in  the  spectral  space  via  a  transfer 
function  as 


Sy{<o)  =  |  H{ o)  1 2Sx(co)  (8) 

where  S,(o))  and  S,(cu)  are  respectively  the  output  and  input  energy-density  spectra  and 
H{(o)  is  the  transfer  function.  The  constant  parameter  linear  system  preserves  the  input 
frequency  at  the  output  and  modifies  the  amplitude  and  phase  of  each  component  of  the 
input  independently.  Similar  types  of  relations  can  be  deduced  for  the  cross-spectral 
quantities.  For  the  case  of  the  open  ocean  wave  field  and  considering  the  pressure  head 
as  input  and  sea  surface  elevation  spectra  as  output,  the  applicable  transfer  function  for 
the  pressure  is 


//»- 


cosh  kz 
cosh  k{h  +  2 ) 


(9) 


which  is  the  same  term  as  in  the  brackets  of  Eq.(7). 

It  was  noted  in  the  formulation  of  the  linear  wave  theory  that  specific  phase 
relations  exist  between  the  sea  surface  elevation  and  pressure.  The  relative  phase  be¬ 
tween  two  quantities  in  spectral  terms  is  defined  as 

(2xy(Oj) 

0°) 


where  Q,,(oj)  is  the  quadrature  spectrum  and  C„,(<y)  is  the  co-spectrum  of  the  two 
quantities. 

The  spectral  transfer  function  relating  pressure  head  to  surface  elevation  as  a 
function  of  frequency  and  depth,  Eq.(9),  can  be  modified  to  include  differential  distance 
and  elevation  between  two  sensors  of  an  array  to  compute  energy  and  cross  -  spectra 
of  surface  slope  in  the  x  -  direction  as 


(11) 

(12) 
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where  A P,  is  subsurface  pressure  difference  between  two  sensors  on  the  x-axis  and  Ajc  is 
differential  distance  along  the  same  axis.  Similarly,  pressure  power  spectra  and  cross  - 
spectra  can  be  computed  along  y-axis. 

B.  THE  SPECIFICATION  OF  A  WAVE  FIELD 

The  complete  description  of  the  wave  field  requires  knowledge  of  frequency  and  di¬ 
rectional  wave  spectra.  Directional  spectra  are  commonly  measured  by  pitch  and  roll 
buoy  and  multi-element  arrays,  which  represent  spectra  measured  at  a  point.  The  meth¬ 
ods  for  computation  of  the  directional  wave  spectrum  are  generally  classified  as 

•  Model  fitting  methods  -  based  on  the  parametric  representation  of  the  spectrum 
as  the  classical  Longuet-Higgins  et  al  (1963)  method. 

•  Model  independent  methods!  Davis  and  Regier,1977),  that  can  be  divided  into  a 
priori  and  data  adaptative  methods. 

A  priori  methods  require  in  advance  a  definition  of  criteria  to  be  verified  without  refer¬ 
ence  to  the  data  in  use.  The  data  adaptative  methods  are  based  in  a  posteriori  assump¬ 
tions  about  the  spectrum  ,  which  depart  from  characteristics  of  the  input  data.  Any  of 
these  processes  of  estimation  are  conditioned  by  the  traditional  constraints  inherent  in 
time  series  analysis  of  statistical  reliability  and  resolution. 

1.  Directional  Spectrum 

The  wavenumber-frequency  characteristic  of  sea  waves  is  contained  in  the  wave 
spectrum  as  given  by  the  Fourier  transform  of  the  correlation  function  R(r ,  t) 

S{K,  <o)  =  — l-rr  R(7,  T)e~i(~K  ’ '  ~a>r)cfr  dr  (13) 

(2rr)  JJ 

This  quantity  is  simply  the  three  dimensional  distribution  of  variance  generally  known 
as  the  power  spectrum.  The  knowledge  of  the  wave  field  at  fixed  positions  (Munk  et  al., 
1963),  permits  the  computation  of  the  correlation  function.  The  definition  of  the  power 
spectrum  as  given  by  Eq.(13)  leads  by  inverse  transforming  to  the  correlation  function, 
that  is 


R(7,  t)  =  J  J S{K,  co)el{K  *'  ~arc)dKd(o 


(14) 
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which  for  the  assumed  stationary  wave  field  is  an  even  function.  Consequently,  the 
power  spectrum  is  real  and  symmetric  about  zero  frequency.  The  correlation-i&mction 
can  then  be  rewritten  as  a  cosine  transform  of  the  power  spectrum 


R(r,  t)  =  J  J S(K,  co)  cos (K,  co)dKdco 


(15) 


The  correlation  function  can  now  be  redefined  as 


R(r ,  r)  =  J  £u)(  cos(ojt)  cos(AT  •  r )  +  sin(a>T)  sin {K  •  7))dKdco  (16) 


The  power  spectrum  can  then  be  calculated  in  terms  of  cosine  and  sine  trans¬ 
forms  of  the  correlation 


S{K,  to)  -  J[C(r ,  to)  +  iQ(7,  co)3e~IK'7d7  (17) 

where  C(r ,  <u)  and  Q(r,  to)  are  respectively  the  co-  and  quadrature  spectra  of  the  corre¬ 
lation.  Such  functions  could  be  calculated  if  the  correlation  is  known  as  continuous 
function  of  the  spatial  coordinate,  although  real  world  observations  lead  to  the  col¬ 
lection  of  data  in  a  small  number  of  positions.  Approximate  methods  of  computing  di¬ 
rectional  spectrum  from  observations  of  the  ,v  ive  field  at  discrete  positions  must  then 
be  considered. 

2.  Classical  Approach 

Lonquet-Higgins,  Cartwright  and  Smith  (1963)  following  the  suggestions  of 
Barber  (1946),  developed  a  method  of  computing  the  directional  spectrum  from  the 
motions  of  a  floating  buoy.  The  method  is  based  on  relationships  between  the  cross 
spectral  quantities  of  the  heave,  pitch,  and  roll  of  the  buoy  and  the  Fourier  coefficients 
of  the  directional  spectrum.  The  directional  spectrum  can  be  represented  as  a  Fourier 
series 


oo 

S{f,9) » -y*  +  cos  nd  +  b„  sin  nd 

nm  1 


(18) 
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where  6  is  the  wave  direction  and  the  coefficients  an  and  bn  are  only  frequency  dependent 
and  represent  the  frequency  distribution  of  energy.  Separating  the  frequency  and  direc¬ 
tional  contributions,  as  it  is  usually  done,  the  directional  spectrum  can  be  rewritten  as 

S(f,6)  =  S(f)D(f,e)  (19) 

where  D(f,$)  is  the  spreading  function.  D(f,6)  corresponds  to  a  unity  area  weighting 
function  that  redistributes  the  energy  contained  in  the  power  spectrum  over  direction, 
ideally  reproducing  real  world  conditions.  The  three  measured  quantities  obtained  with 
a  buoy  can  be  directly  related  with  the  harmonic  representation  of  the  sea  surface.  The 
heave  corresponds  to  the  instantaneous  sea  surface  elevation  and  the  pitch  and  roll  are 
associated  with  the  spatial  partial  derivatives  of  the  surface  elevation.  Considering  a 
single  harmonic  of  the  sea  surface  representation 

ri(t)  =  a  cos{K  •  x  —cot)  (20) 

the  two  horizontal  cartesian  components  of  the  sea  slope  can  be  written  as 


cri(t) 

Vx(t)  = =  {  ~iK  cos  6}q(t) 

(21) 

dri(t) 

Vy(')  " - {-iK  sin  d)r,(t) 

(23) 

Cross  spectral  quantities  can  be  computed  for  the  instantaneous  surface  ele¬ 
vation  and  surface  slope,  and  related  with  the  Fourier  coefficients  of  the  directional 
spectrum.  These  cross  spectral  quantities  can  be  represented  as 


C^if)  =  Si/)!2*  Dtf,d)dd 
Jo 


(23) 


Q^,jj)=KS{f)  r 

Jo 

Q  <J)  =  KS(f)  r 
Jo 


D{f,9)  cos  Odd 


D{f,6)  sin  Odd 


(24) 


(25) 
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(26) 


C„  „,(/)  =  K2S(J)  I* Dlffi)  cos 2ddd 
j0 

C,  ,  (/)  -  tf2S(/)  f  UD(f,6)  sin26dd  (27) 

C  (/)  =  tf2S(/)  f  2,rZ)(/:0)  cos  6  sin  M0  (28) 

*  ’  Jq 

Considering  the  directional  spectrum  in  terms  of  its  first  five  Fourier  coefficients,  that  ic 
S(f,6)  =  +  a,  cos  -  6,  sin  6  +  a2  cos  26  +  b2  sin  26  (2‘ri 

,  it  can  be  shown  that  these  coefficients  are  given  in  terms  of  the  cross  spectral  quaiu 
ties,  Eqs.(23-28),  by 


„  \  r 

°o  —  n  t'n*i 

(30) 

=  Tk  Q’”'* 

(31) 

(32) 

a2  ~  ^^2 

(33) 

b2  =  — C„  „ 
nK 2 

(34) 

and  equivalently,  if  a  normalization  by  auto-spectrum  is  done,  the  Fourier  coefficients 
of  the  unit  area  spreading  function  are  determined.  This  approach  constitutes  the  clas¬ 
sical  LCS  method.  The  approximate  representation  of  the  directional  spectrum  in  terms 

of  its  Fourier  series  truncated  to  its  first  five  coefficients  leads  to  a  directional  distrib- 

ution  of  energy  as  given  by  a  cosine-bell  function.  Such  a  representation  is  of  limited 
resolution,  and  typically  has  large  negative  side  lobes  (Figure  2V 
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To  solve  the  problem  of  the  unrealistic  negative  side  lobes,  LCS  imposed  a 
constraint  of  positivity  on  the  directional  spreading  function  which  resulted  in  a 
smoothing  effect.  However  this  approach  has  the  consequence  that  the  resolution  of  a 
distribution  narrower  than  cos4(  — )  is  not  possible.  It  must  be  noted  that  this  artifice 
does  not  solve  the  real  problem  that  is  related  to  unresolved  harmonics  of  the  severely 
truncated  Fourier  series  representation.  Also,  the  LCS  approach  is  not  capable  of  re¬ 
presenting  multi-modal  distributions  of  energy  that  are  frequently  observed  in  the  real 
ocean,  as  can  be  concluded  from  the  truncated  Fourier  series  given  by  Eq.(29)  and 
shown  in  ( Figure  2). 

The  method  just  described  can  also  be  applied  to  measurements  using  slope  ar¬ 
rays,  w’hich  consist  of  multi-point  arrangements  of  pressure  transducers.  Typically  the 
sensors  in  the  slope  arrays  are  arranged  in  a  square.  The  real  slope  is  approximated  by 
dividing  the  differential  pressure  of  a  pair  of  sensors  by  their  separation  while  assuming 
the  sea  surface  slope  as  constant  between  each  pair  of  sensors.  The  method  of  com¬ 
puting  the  directional  distribution  of  energy  then  proceeds  as  before  for  the  case  of 
measurements  with  a  floating  buoy.  The  sea  surface  elevation  is  obtained  from  the  time 
series  of  pressure  by  applying  the  appropriate  transfer  function.  The  approximation  of 
the  surface  slope  introduces  errors  in  the  computation  given  by  (Seymour  and  Higgins, 
1977) 


dr\  sinffffeoq)  dr] 

‘  =  Kb*, 

2 


(35) 


where  a,  =  (  cos  6,  sin  6)  for  i=  1,2  and  b  is  the  horizontal  sensor  separation.  It  can  be 
seen  that  the  error  of  such  an  approximation  is  a  function  of  the  wavenumber  vector  and 
the  size  of  the  array,  decreasing  in  magnitude  as  these  quantities  are  reduced.  It  is  then 
natural  to  employ  arrays  as  small  as  possible,  limited  by  the  resolution  and  accuracy  of 
the  pressure  sensors. 

3.  Exact  Fourier  Coefficient  Representation  Method 

The  LCS  approach  provides  a  way  of  computing  the  first  five  Fourier  coeffi¬ 
cients  of  the  series  expansion  of  the  directional  spreading  function.  These  coefficients  are 
obtained  via  cross-  spectral  analysis  of  triorthogonal  components  of  the  wave  field 
measurements.  The  number  of  Fourier  coefficents  are  limited  because  of  the  measure¬ 
ment  technique  normally  used,  floating  buoy  or  slope  array  of  four  sensors  in  a  square. 
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Figure  2.  Typical  spreading  function  of  the  LCS  method:  Small  directional  resol¬ 
ution  and  unrealistic  negative  side  lobes. 

Such  a  system  can  only  give  information  of  the  wave  field  through  second  order  terms. 
The  coelficients  determined  can  be  represented  by  an  infinite  number  of  dilferent  dis¬ 
tributions  whose  Fourier  expansion  is  common.  Examples  of  distributions  having  the 
same  ouricr  coefficients  can  be  seen  in  Figure  3.  It  is  then  natural  to  choose  a  phys- 
icall;  calisable  distribution,  which  matches  the  coefficients  determined  by  _ observation 
and  more  exactly  simulates  the  true  directional  distribution. 

Grauzini$(1989)  exploits  the  matching  of  a  set  of  Fourier  coefficients  by  dif¬ 
ferent  distributions  to  develop  a  new  method  of  computing  directional  spectra.  The 
previously  defined  directional  spreading  function  D(0)  can  be  represented  as 


m  -  ( 1  +  cost*#,)  +  b„  sinfM*)]}  (37) 

where  ak  and  bt  are  the  unitary  Fourier  coefficients.  Such  coefficients  enter  the  Tinition 
of  the  polar  coefficients  c*  and  6t  as  given  by 
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(37) 


Ck  =  ak  +  ibk  =  ckelk9k 

Any  symmetric  unity  area  function  can  be  determined  solely  by  its  canonical  coefficients 
defined  as 


mk  =  X(6)  cos  kddd  (38) 

where  X(6)  is  any  arbitray  unit  area  function  symmetric  around  0  =  0.  This  represen¬ 
tation  is  simply  the  cosine  transform  of  an  even  function.  It  is  seen  that  a  relation  can 
be  established  between  the  canonical  coefficients  and  the  Fourier  coefficients. 

ak  +ibk  =  ckelk8>  =  mkelkpk  (39) 

where  the  polar  form  is  now  used  for  the  Fourier  coefficients  of  the  spreading  function. 
Using  linear  superposition,  any  unit  area  function  can  be  reproduced  by  a  weighted  sum 
of  symmetric  unit  area  density  functions  oriented  at  different  directions  as  given  by 

ck*lMi  =  YjWJmk^e‘kPl  (4°) 

J 

where  w ,  is  the  weighting  factor  and  /?y  is  the  orientation  of  each  of  the  components 
considered  in  the  series.  Eq.(40)  constitutes  the  generating  equation  for  the  polar  coef¬ 
ficients  where  on  the  right  hand  side  wJt  m*(/)  and  /?,  are  unknowns.  The  left  hand  side 
of  Eq.(40)  is  determined  from  the  cross  spectral  quantities  of  the  wave  field  typically 
known  through  second  order.  A  system  of  equations  can  be  established  as 


mx[w^eip'  +w2e^J]  =  cxeWl 

(41) 

+w2en^J]  =  c2e^J 

(42) 

which  can  be  seen  to  be  underdetermined.  It  must  be  noted  that  in  the  system  Eq.(41-42) 
it  is  assumed  that  the  canonical  coefficients  are  common  to  both  of  the  distributions  of 
energy  considered.  This  was  the  approach  followed  by  Grauzinis  (1989)  to  obtain 
analytic  solutions  for  the  system  of  generation  equations.  By  the  condition  of 
orthogonality  of  the  exponentials  of  different  order  and  the  side  condition  that  the  sum 
of  the  weights  must  be  equal  to  unity,  analytic  solutions  are  formed.  The  analytic 
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DIRECTION 


I  iguic  3.  Distribution  of  energy:  The  three  distribution  have  a  common  Fourier 
expansion. 

solutions  constitute  an  implicit  underdetermined  system  of  four  equations.  The 
unknowns  of  this  system  are  the  two  canonical  coefficients,  the  orientation  of  the  two 
unit  area  density  functions  and  one  of  the  weighting  functions.  To  solve  this  system, 
further  constraints  are  needed.  Assuming  a  functional  relationship  between  the 
canonical  coefiiceints  m,  and  m,  particular  solutions  are  obtained.  The  original  relation¬ 
ship  adopted  was  a  power  law  between  the  two  canonical  coefficients 

m2  =*  m\  (43) 

that  for  different  values  of  r  originated  different  distributions  of  energy  as  analytic  sol¬ 
utions.  Such  distributions  are  of  decreasing  sharpness  as  the  magnitude  of  the  exponent 
increases.  For  the  values  of  the  exponent  2,  3  and  4,  the  resulting  distributions  in  its 
unimodal  form  are  represented  in  Figure  4  using  two  different  orientations  and  weights 
also  allows  this  method  to  represent  bimodal  distributions  of  energy.  The  variable 
sharpness  of  the  model  distributions  considered  will  permit  the  matching  of  observed 
distributions  of  considerably  more  peakedness  than  the  classical  LCS  method.  For  a 
canonical  coefficient  of  .9,  the  haft-power  width  of  the  distribution  D4  is  12  degrees,  that 
is  about  one-third  of  the  typical  resolution  of  the  cosine-bell  distribution. 
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Knowledge  of  Fourier  coefficients  of  order  higher  than  the  second  introduces 
considerable  advantage  to  this  method.  Some  of  the  constraints  can  be  dropped,  such 
as  the  functional  relationship  of  the  canonical  coefficients.  Moreover,  the  matching  of 
coefficients  of  higher  order  increases  the  resolution  of  the  method. 

The  exact  matching  of  the  Fourier  coefficients  through  second  order  of  any 
spreading  function  can  be  obtained  by  the  method  just  described.  The  problem  of  what 
criteria  to  use  for  the  selection  of  the  distribution  that  best  matches  the  real  data  must 
also  be  considered.  It  should  be  emphasized  that  there  exist  an  infinity  of  distributions 
that  can  be  considered  to  match  the  same  set  of  Fourier  coefficients.  The  problem  of 
selecting  the  best  fitting  distribution  partly  amounts  to  the  matching  of  the  ocean 
peakedness.  Cartwright  (1963)  considered  this  problem  when  computing  directional 
wave  spectra  by  the  LCS  method.  He  concluded  that  information  could  be  obtained 
from  the  ratio  of  the  magnitude  of  the  coefficients  C„,  Eq.(37).  The  magnitude  of  such 
coefficient  is  given  by 

I  Q  I  =  +  b\  (44) 

Recalling  the  functional  relationship  of  the  canonical  coefficients  introduced  to  allow  for 
the  solution  of  the  generating  system,  Eq.(43),  it  is  concluded  that  in  the  Grauzinis 
method  it  is  more  natural  to  consider  the  ratio  of  the  natural  logarithms  of  the  C*  coef¬ 
ficients,  that  is 


LogRa,io~HQ)  ‘4S> 

which  represents  the  greatest  power-law  of  the  canonical  coefficients  that  can  match  the 
input  Fourier  coefficients.  Cartwright  (1963)  verified  that  the  ratios  obtained  from  the 
observed  data  were  consistently  greater  than  the  values  corresponding  to  the  cosine-bell 
function  of  the  LCS  method.  It  was  hypothesized  that  such  difference  could  be  due  to 
the  multimodality  of  the  wind- wave  spectrum.  Recalling  that  the  ratio  of  the  coefficients 
C*  as  an  indicator  was  deduced  from  uni-modal  considerations;  it  does  not  seem  correct 
to  make  use  of  it  to  infer  of  the  adequacy  of  multi-modal  distributions.  Grauzinis  (1989) 
points  out  that  the  natural  logarithimic  ratio  will  exceed  the  power  law  of  canonical  co¬ 
efficients  increasingly  as  the  beam  separation  gets  larger.  Moreover,  noise  affects  the 
coefficients  C2  and  C,  differently.  Eqs.(23-28)  show  that  the  effects  of  noise  are  cancelled 
for  the  coefficient  Oj ,  and  consequently  the  Log  Ratio  will  be  decreased. 
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Figure  4.  Three  model  distribution:  Uni-modal  representation  of  three  different 
model  distributions  for  two  values  of  the  canonical  coelficient.  D2  ,D3 
ind  D  orrespond  respectively  to  the  exponent  2,3  and  4  for  the  func¬ 
tional  i  tionship  of  the  canonical  coefficients  Eq.(43).  The  sharpness 
of  the  distribution  decreases  with  the  magnitude  of  exponent.  For  a 
value  of  the  canonical  coefficient  equal  to  .9  the  half  power  width  of  D4 
is  12  degrees. 

C.  WAVE  SPECTRUM  TRANSFORMATION 

According  to  Mehaute  and  Wang  (1982),  all  the  formulations  of  wave  spectrum 
transformation  over  a  shoaling  bottom  are  based  on  the  premise  that  the  wave  energy 
associated  with  a  narrow  frequency  band  stays  within  that  band  during  the  transforma¬ 
tion  and  that  it  obeys  linear  superposition.  (The  shift  of  energy  from  frequency  to  fre¬ 
quency  is  considered  to  be  the  results  of  nonlinearity.)  For  each  frequency,  the  energy 
level,  which  is  characterized  by  the  square  of  a  free  surface  elevation,  >?J(f),  is  treated  as 
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an  invariant  during  the  transformation.  Therefore,  the  law  of  tansformation  which  has 
been  established  for  monochromatic  waves  is  applied  as  if  each  value,  tj( 0,  were  identical 
to  a  periodic  wave  of  the  same  amplitude  and  frequency,  and  that  the  energy  contained 
in  each  band  and  direction  travels  along  its  corresponding  wave  ray  at  group  velocity. 

The  first  rigorous  theoretical  approach  to  the  evolution  of  the  wave  spectrum  is  due 
to  Longuet-Higgins  (1957),  who  demonstrated  that  the  energy  spectral  density  in  the 
wave  number  space  remains  constant  as  one  follows  any  wave  group  along  a  wave  ray 
(Phillips  1966).  The  demonstration  is  based  on  the  conservation  of  energy  and  the  law 
of  wave  refraction.  It  is  shown  that  S(kltk2)  =  constant  along  a  wave  orthogonal,  or 
~  S(ku  k2)  =  0  .  Other  demonstrations  have  also  been  presented  by  Karlsson  (1969), 
Collins  (1972),  and  Krasitskiy  (1974)  in  analogy  with  geometrical  optics  and  the  ray 
particle  analogy. 

Following  their  approach,  let  us  consider  a  wave  spectrum  defined  in  terms  of  space 

(jc,>')<  time(r),  and  wave  number  (kltk2)  as  S(klt  k2,  x,y,  t).  By  virtue  of  the  aforemen- 
dS 

tioned  principle  (-^—  =  0),  then 

cS  .  cS  dx  ,  cS  dy  .  8S  dk\  .  8S  dki  „ 

— - h  -z - —  +  — - —  +  — - - - h  -rr - T“  =  0  (46) 

cl  dx  dl  cy  dl  ckx  dt  ok2  di 

The  first  term  is  generally  taken  equal  to  zero  under  steady-  state  condition.  The  second 
and  third  terms  are  the  convective  terms  and  are  due  to  the  variation  of  energy  level  as 
a  result  of  the  spatial  variation.  The  two  last  terms  are  due  to  shoaling  and  refraction, 
is  the  velocity  of  energy  propagation  in  the  x  direction  and 


(47) 

Similarly: 

±=j2._Vsiae-vy  (48, 

where  V:  group  velocity,  k:  wave  number.  Since  w2  =  gk  tanh  kh,  a>  is  a  function  of  k 
and  h  .  The  depth,  h,  depends  exclusively  on  (jt,^),  thus  we  can  write 
to  =  co[k,/i,(jc,y)]  =  (o(kj,y).  Then 


do  _  do  dk\  do  dx  do  dk2  do  dy 
dt  ~  ckx  dt  dx  dt  ck2  dt  dy  dt 
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in  which  C  *  the  phase  velocity,  C  «•  — It  is  recalled  that  the  ray  equations  are  (Munk 
and  Arthur  1951) 


(56) 


de 

ds 


dx  0  dy  . 
j  =  cos  — r-  —  sin  0 
ds  ds 


1 


•  f  sin  6  -  cos  d 

\  ex  cy 


in  which  s  =  the  distance  along  a  ray.  Inserting  these  into  Eq.(55)  yields  simply: 

■^fCCKSWO-O  (57) 

i.e.: 

CVS(f,6)  —  constant  (58) 


along  a  wave  ray.  Note: 


dS(kuk2) 

dt 


_V _ d_ 

Into  ds 


[CFS(/;0)]=O 


(59) 


The  result  obtained  in  Eq.(59)  is  the  same  as  Longuet-Higgins.  Eq.(58)  can  be  written 
further  as 


—  S{ffirX,y)  —  constant 


or 


S(W,fi)=^S0(W,90)  (61) 

where  the  subscript  refers  to  deep  water  quantities.  In  other  words,  the  ratio  between 
shallow  water  DSD  and  deep  water  DSD  is  a  universal  function  of  depth  and 
wavenumber  where 

(62) 

Figure  5  shows  the  plots  and  numerical  values  of  the  ratio  between  deep  and  shallow 
DSD  Eq.(62)  at  Marina  and  Santa  Cruz  versus  different  frequencies. 
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Figure  5.  Transformation  constant  along  a  ray  path 
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III.  DATA  ACQUISITION  AND  SENSORS 


A.  LOCATION 

Monterey  Bay  was  selected  as  the  study  area,  which  poses  complex  problems  due  to 
the  highly  irregular  bathymetry  of  the  Monterey  submarine  canyon.  A  deepwater  wave 
buoy  and  two  shallow  water  slope  arrays  were  utilized  for  wave  data  acquisition.  The 
locations  of  sensors  are  listed  in  Table  1,  see  Figure  1  for  bathymetry  and  location. 


Table  1.  LOCATIONS  OF  SENSORS 


No. 

Type  of  Sensor 

Location 

Depth  (M) 

1. 

NDBC  roll/ pitch  wave 
buoy 

36°48.0'N, 

122°23.0'59.0"W 

2003 

2. 

Array  of  pressure  gauges 
(4)  Marina 

36°42.0'N,  121°48.9'W 

15.0 

3. 

Array  of  pressure  gauges 
(4)  Santa  Cruz 

36°57.0'N,  122°0.2'W 

13.1 

In  comparing  coincident  data  sets,  on  the  average,  a  time  lag  of  one  hour  was  assumed 
for  the  waves  to  travel  from  the  deep  water  buoy  to  the  shallow  water  array  sites. 

B.  DEEP  WATER  WAVE  DATA 

Deep  water  wave  data  of  the  pitch/roll  buoy  were  obtained  from  the  NDBC  (Na¬ 
tional  Data  Buoy  Center).  The  NDBC  3-meter  buoy  measures  heave,  pitch,  and  roll 
movements  with  respect  to  a  3-axes  magnetometer.  The  magnetometer  axes  are  aligned 
with  fore  and  after  (bow  positive),  athwartship  (starboard  positive)  and  vertical  (positive 
upward)  axes  of  a  buoy,  see  Figure  6. 

The  analog  sensor  measures  vertical  acceleration  and  displacement,  hull  pitch  and 
roll.  The  magnetometer  bow  and  starboard  output,  together  with  vertical  displacement, 
pitch  and  roll  output  from  analog  sensors  are  sent  to  the  Directional  Wave  Data  Ana¬ 
lyzer  (DWDA)  after  voltage  amplification  and  filtering. 

The  hull  azimuth  angle  is  calculated  from  bow  and  starboard  magnetic  components 
with  necessary  corrections  made  for  pitch,  roll,  and  hull  magnetic  fields.  The  magnetic 
azimuth  is  converted  into  true  azimuth  using  the  magnetic  variation  for  the  site.  The 
true  azimuth  is  then  used  to  calculate  East-West  and  North-South  components  of  hull 
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Figure  6.  Floating  Buoy 


slope.  One  hertz  time  record  of  hull  displacement  and  slope  are  stored  in  DWDA  for 
spectral  analysis  and  transmission  to  shore  station  through  GOES  (Steele  et  al.,  1985). 

Data  provided  on  magnetic  tape  included  environmental  parameters,  wave  spectra, 
directional  wave  parameters,  and  co  and  quad-spectra.  Directional  wave  parameter  data 
were  utilized  to  determine  directional  wave  spectrum. 

C.  SHALLOW  WATER  WAVE  DATA 

The  calibrated  data  for  the  presst  2  gauge  arrays  were  obtained  from  the  Scripps 
Institute  of  Oceanography.  Array  consists  of  four  pressure  gauges  located  in  a  six-meter 
square  area  which  measures  pressure  head  in  centimeters  of  water,  See  Figure  7.  Each 
sensor  is  located  about  50  to  100  cm  above  the  bottom.  Pressure  power  spectrum  is 
calculated  which  in  turn  is  transformed  into  surface  elevation  rj  and  surface  slopes  tj,  and 
rif  spectra  by  applying  linear  wave  theory  transfer  function  as  described  in  Chapter  II. 
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Figure  7.  Slope  Pressure  Array 


For  the  computation  of  the  cross-spectra,  the  following  values  were  utilized: 

•  sampling  interval  =  At  *»  2  sec 

•  total  data  points  «  2048 

•  data  points  for  one  sample  *  64 

•  nyquist  frequency  *  0.25  Hz 

•  frequency  bandwidth  *  bf  =*  0.0078125  Hz 

The  spectra  were  cut  at  a  frequency  of  0.2  Hz  in  the  higher  frequency  range  due  to 
overestimation  by  the  transfer  function.  The  wave 'period  for  high  quality  estimate  is 
approximately  5-20  seconds. 


D  WAVE  DATA  SELECTION 

Selected  wave  cases  from  January  19S8  were  used  in  the  analysis.  During  the  winter 
season,  weather  conditions  were  generally  favorable  for  creating  higher  energy  waves. 
The  selected  days  are  3  January  0800  PST,  18  January-  at  0200  PST  and  29  January-  at 
2100  PST.  These  periods  provide  two  cases  of  narrow  band  energy  spectra  at  different 
frequencies  and  directions  and  a  third  case  of  a  major  storm  event.  The  California  coast 
was  hit  by  an  unusually  severe  storm  during  16-18  January  1988  which  resulted  in  high 
waves,  strong  winds,  and  sea  level  well  in  excess  of  predicted  values  (Cavan,  et  al.,  1988). 
Widespread  damage  occured  along  the  southern  California  coast.  This  storm  developed 
about  300  \M  west  of  San  Francisco  on  16  January  1988  traveled  south-east  at  about 
33  knots.  The  storm  center  passed  off  Monterey  Bay  on  the  morning  of  17  January  1988 
continuing  on  its  SE  track  and  making  landfall  at  Avila  Beach  in  Central  California. 
The  energy  spectra  for  these  cases  are  depicted  in  Figure  8,  9,  and  10  comparing  the  deep 
water  buoy  and  shallow  water  pressure  gauge  measurements  at  Marina  and  Santa  Cruz. 
Deep  water  wave  characterics  of  selected  cases  are  given  in  Table  2. 


Table  2.  DEEP  WATER  WAVE  CHARACTERISTICS  OF  SELECTED  CASES 


Characteristics 

3  Jan  88 

18  Jan  88 

29  Jan  88 

tf,,. 

1.87  m 

5.12  m 

2.65  m 

0.13  Hz 

0.07  Hz 

0.09  Hz 

EBI^H 

0.18  Hz 

0.10  Hz 

0.12  Hz 

Direction „„„ 

o 

O 

oo 

C4 

300° 

275° 

Reviewing  Table  2  and  other  published  data,  deep  water  waves  generally  travel 
eastward  during  the  month  of  January  with  slight  variability  of  about  20°  on  eigher  side 
of  true  East.  This  direction  is  dominant  in  the  energetic  middle  frequency  range  of  0.08 
to  0. 1 3  Hz.  At  other  frequencies,  the  direction  tends  to  be  variable. 

The  significant  wave  heights,  calculated  as  H,  —  4yja2,  where  a2  is  variance  calcu¬ 
lated  as  the  area  under  the  energy  density  spectra,  are  given  in  Table  3. 

E.  TRANSFORMATION  OF  DEEP  WAVE  SPECTRUM 

The  directional  spectrum  density  at  wavenumber  space  is  constant  along  ray  path. 
Transforming  from  wavenumber  space  to  frequency  and  direction  space,  the  shallow 
water  DSD  is  related  to  the  deep  water  DSD  by  a  universal  function  of  depth,  h,  and 
wavenumber,  k,  or  frequency,/,  Eq.(62),  however,  the  relationships  between  the  shallow 
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S(f),  M**2/llz 


Frequency,  Hz 


Figure  8.  Measured  spectral  densities  on  3  January  1988,  0800  PST 
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Figure  10.  Measured  spectral  density  on  29  January  1988,  2100  PST 
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Table  3.  SIGNIFICANT  WAVE  HEIGHT  (UNIT:  METER) 


Buoy 

Marina 

Santa  Cruz 

3  Jan  1988 

1.87  M 

0.8  M 

0.9  M 

18  Jan  1988 

5.12  M 

5.15  M 

2.4  M 

29  Jan  1988 

2.65  M 

1.77  M 

0.94  M 

water  angle  and  deep  water  angle  are  still  to  be  found.  The  Dobson  linear  refractir 
model  was  used  for  this  purpose.  The  reverse  projection  method  was  employed.  That  is, 
the  linear  refraction  model  program  was  run  backward  from  shallow'  water  location  to 
deep  water,  beginning  at  the  shallow  water  locations  at  Marina  and  Santa  Cruz.  Figure 
12  -  15.  show  rays  at  the  particular  frequencies  of  0.03,  0.06,  0  ).  0.13  Hz  from  .arina 
and  Santa  Cruz  propagated  offshore  in  increments  of  0.5  degrees  over  the  range  of  all 
possible  incoming  wave  angles.  Once  the  de.  p  water  is  reached,  the  rays  are  stopped 
and  deep  water  wave  directions  are  recorded.  As  can  be  seen  the  low  frequency  waves 
are  more  refracted  by  the  bathymetry  than  high  frequency  w'aves. 

An  example  of  the  shallow  water  angle  versus  deep  water  angle  at  Marina  for  0.06 
Hz  is  show-n  in  Figure  15.  The  doited  line  is  the  least  square  first  order  fit,  where  a  is 
the  slope  and  b  is  the  intersection  at  shallow  water  angle  axis  when  deep  water  angle 
becomes  zero. 

The  slope  and  intercept  of  least  squa.  inear  fits  at  various  frequencies  and  lo¬ 
cations  are  given  in  Table  4  and  plots  given  *.  Appendix  A.  As  frequency  increases,  the 
slope  increases  to  one  and  the  offset  decreases  toward  zero.  Since  slope  less  than  one 
amounts  to  shrinking  the  deep  water  angular  range,  it  indicates  that  smaller  slope  will 
concentrate  wave  energy  into  a  smaller  direction  range  and  decrease  the  spectral  density. 
The  offset  is  equivalent  to  rotating  the  deep  water  DSD  to  a  different  direction. 
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Figure  12.  Ray  traces  from  Marina  to  deep  water 
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X 

Low  Frequency 

Figure  13.  Ray  traces  from  Santa  Cruz 

4 


3 


Shallow  Water  Angle  vs  Deep  Water  Angle  at  Marina,  0.06  Hz 


Figure  15.  Example  of  the  shallow  water  angle  versus  deep  water  angle 


Table  4.  EFFECTS  OF  LINEAR  REFRACTION  _ 

Marina  Santa  Cruz 


Frequency 

slope 

offset 

slope 

offset 

0.03  Hz 

0.1164 

242.0 

0.3290 

126.9 

0.06  Hz 

0.3020 

1S9.0 

0.3974 

116.6 

0.09  Hz 

0.4165 

157.9 

0.4998 

98.0 

0.13  Hz 

0.6938 

82.8 

0.6824 

62.1 

IV.  RESULTS  AND  DISCUSSION 


A.  COMPARISON  OF  LCS  AND  EFC  MEASURED  DIRECTIONAL  SPECTRAL 
DENSITIES 

A  flow  chart  of  the  calculation  of  measurements  is  shown  in  Figure  16.  The 
TRIDEN  program  calculates  the  EFC  spectrum  and  compares  it  with  the  spectrum 
measured  in  shallow  water.  The  directional  distribution  of  measured  waves  in  polar  plot 
form  are  shown  in  Figures  17, 18,  and  19  at  the  three  locations,  Buoy,  Marina,  and  Santa 
Cruz  on  3  January  1988.  These  axis  are  aligned  with  the  earth  directions.  The  deep 
water  wave  approach  from  northwest  or  west.  Four  frequencies,  0.03, 0.06,  0.09  and  0.13 
Hz  are  examined.  At  Marina,  waves  come  mostly  from  west,  and  at  Santa  Cruz,  from 
southwest.  The  plots  show  the  energy  is  more  widely  distributed  in  deep  water  and  that 
the  spectrum  narrow's  in  shallow  water,  except  in  the  lowest  frequency. 

The  energy  density  distribution  as  a  function  of  frequency  and  direction  are  shown 
in  Figures  20  -  31,  for  3  January  1988  at  deep  water  Buoy,  Marina,  and  Santa  Cruz. 
Both  of  the  LCS  method  and  EFC  method  indicate  the  same  directionality  of  measured 
data,  but  EFC  method  shows  a  much  sharper  peak  than  LCS  method.  Also  LCS  results 
sometimes  show  negative  side  lobes,  which  do  not  exist  in  the  EFC  results. 

The  wave  energy  spectra  at  frequency  0.13  Hz,  are  shown  in  Figures  32,  33,  and  34. 
Again  ,  the  EFC  method  results  in  much  higher  resolution.  For  this  reason,  only  the 
EFC  are  discussed  for  other  results.  The  other  two  days,  18  and  29  January  1988,  3-D 
and  2-D  also  show  similar  results.  All  graphs  are  given  in  Appendix  B.  The  deep  water 
buoy  and  shallow  water  array  (Marina  and  Santa  Cruz)  energies  are  compared  in  Table 
5,  6  for  the  different  methods.  As  can  be  seen  the  LCS  and  EFC  methods  give  the  same 
total  energies  under  the  spectra. 


Table  5.  DEEP  WATER  BUOY  ENERGY  (AP) 


DATE 

EFC 

LCS 

D2 

COMPOSITE 

3  Jan  88,  0800 

0.22 

0.22 

0.22 

18  Jan  88,  0200 

1.642 

1.641 

1.642 

21  Jan  88,  2100 

0.443 

0.442 

0.443 
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Figure  16.  The  flow  chart  of  measurement  calculation 


Table  6.  SHALLOW  WATER  ARRAY  ENERGY  (A P) 


DATA 

EFC 

LCS 

D2 

COMPOSITE 

Marina 

Santa 

Cruz 

Marina 

mmrm 

■33131 

Marina 

Santa 

Cruz 

3  Jan  88.0800 

0.042 

0.054 

0.042 

0.054 

0.042 

0.054 

18  Jan  88. 
0200 

1.659 

0.363 

1.656 

0.362 

1.656 

0.362 

21  Jan  88. 
2100 

0.196 

0.057 

0.196 

0.057 

0.196 

0.057 

B.  REFRACTION  MODEL  PREDICTIONS 

The  deep  water  spectrum  is  transformed  to  shallow  water  using  Eq.(62).  The  shal- 
low  water  spectrum  is  simply  equal  to  deep  water  spectrum  multiplied  shoaling. 
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..  EFC  SPECTRAL  DENSITY  (M*^/HZ/RAQ) 


BUOY 


Figure  17.  The  energy  direction  of  deep  water  on  3  Jan  1988. 


37 


38 


SANTA  CRUZ 


Figure  19.  The  energy  direction  of  shallow  water  at  Santa  Cruz  on  3  Jan  1983. 
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Figure  20.  The  LCS  directional  spectral  density  in  3-D  plot:  On  3  Jan  1988  at 

* 

deep  water  Buoy. 
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directional  spectra 

Deep  Water  Buoy 
Jan  03,  19e8  0800  POT 


Figure  21.  The  EFC  directional  spectral  density  in  3-D  plot:  On  3  Jan  1988  at 
deep  water  Buoy. 
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Figure  22.  The  LCS  directional  spectral  density*  in  contour  plot:  On  3  Jan  1988 

at  deep  water  Buoy. 

* 


Figure  23.  The  EFC  directional  spectral  density  in  contour  plot:  On  3  Jan  1988 
at  deep  water  Buoy. 

e 
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DRECTC'.'AL  SPECTRA 
Menne  -w  Wcter  Amoy 
Jan  C?  0800  POT 


Figure  24.  The  LCS  directional  spectral  density  in  3-D  plot:  On  3  Jan  1988  at 
Marina. 
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Figure  25. 


The  EFC  directional  spectral  density  in  3-D  plot:  On  3  Jan  1988  at 

* 

Marina. 
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Figure  27.  The  EFC  directional  spectral  density  in  contour  plot:  On  3  Jan  1988 
at  Marina. 

a 
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Figure  28.  The  LCS  directional  spectral  density  in  3-D  plot:  On  3  Jan  1988  at 

4 

Santa  Cruz. 


Figure  29. 


The  EFC  directional  spectral  density  in 
Santa  Cruz. 
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Figure  31.  The  EFC  directional  spectral  density  in  contour  plot:  On  3  Jan  1988 
at  Santa  Cruz. 
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Figure  32.  The  LCS  &  EFC  directional  spectral  densities  in  2-D  plot:  On  3  Jan 
1988  at  deep  water  Buoy. 
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Figure  33.  The  LCS  &  EFC  directional  spectral  densities  in  2-D  plot:  On  3  Jan 
1988  at  Marina. 
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Figure  34.  The  LCS  &  EFC  (directional  spectral  densities  in  2-D  plot:  On  3  Jan 
1988  at  Santa  Cruz. 

* 
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Table  7.  COMPARISON  OF  TRANSFORMED  AND  MEASURED  DATA  DIF¬ 


FERENCE  IN  DEGREE 


Location 

Marina 

Santa  Cruz 

Frequency 

0.06 

Hz 

0.09 

Hz 

0.13 

Hz 

0.06 

Hz 

0.09 

Hz 

0.13 

Hz 

03  Jan  19S8 

25° 

0° 

0° 

15° 

0° 

30° 

18  Jan  1988 

30° 

95° 

50° 

5° 

15° 

25° 

29  Jan  1988 

5° 

0° 

5° 

40° 

5° 

5° 

Table  8.  MEASURED  AND  TRANFORMED  ENERGY  SPECTRUM  DATA  AT 


MARINA  (API HZ) 


FREQ. 

0.06  Hz 

0.09  Hz 

0.13  Hz  | 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

3  Jan  88.  0800 

0.026 

0.123 

0.062 

0.078 

0.401 

0.449 

18  Jan  88.  0200 

12.8 

3.85 

1.54 

2.46 

4.11 

2.90 

29  Jan  88.  2100 

0.103 

0.103 

1.28 

1.28 

0.180 

0.205 

Table  9.  MEASURED  AND  TRANFORMED  ENERGY  SPECTRUM  DATA  AT 
_ SANTA  CRUZ  ( API  HZ) 


FREQ. 

0.06  Hz 

0.09  Hz 

0.13  Hz  | 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

MEAS¬ 

URED 

TRANS¬ 

FORMED 

3  Jan  88.  0800 

0.072 

0.187 

0.045 

0.118 

0.205 

0.277 

18  Jan  88.  0200 

2.87 

3.25 

1.03 

2.00 

0.411 

1.54 

29  Jan  88.  2100 

0.062 

0.16 

0.616 

1.03 

0.168 

0.205 

V  v 

(-y-),  and  dispersion,  (-g-)  and  the  energy  density  is  replotted  based  on  the  angle  con¬ 
version  between  deep  and  shallow  water,  determined  by  the  linear  wave  refraction  model. 

The  Marina  location  is  shadowed  by  Point  Pinfis  in  the  southwest  and  Point  Santa 
Cruz  in  the  northwest.  High  frequency  waves,  O.ll  Hz  and  above,  are  limited  to  di¬ 
rections  from  235°-310°  at  the  array  site,  whereas  low  frequency  waves  may  arrive  from 
even  southerly  directions  after  considerable  refraction  (refer  to  Figures  11,  12). 
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The  Santa  Cruz  site  poses  more  complex  problems  than  Marina  because  of  its  lo¬ 
cation.  The  waves  approaching  from  the  westerly  direction  -  ke  almost  a  90®  turn  to 
arrive  at  the  shoreline.  The  low  frequency  waves  travelling  :  rth  and  north-eastward 
may  arrive  at  this  location.  The  ray  patterns  show  (Figures  13,  14)  scattered  and  cross¬ 
ing  rays  (caustics)  from  the  west  or  higher  approach  angle  at  low  frequency.  The  ray 
patterns  become  more  regular  with  increasing  frequency. 

C.  COMPARISON  OF  DEEP,  SHALLOW  AND  TRANSFORMED 
DIRECTIONAL  SPECTRAL  DENSITIES 

Transformed  and  measured  shallow  wave  directional  spectra  for  Marina  are  com¬ 
pared  in  Figure  35  -  37.  On  3  January  1988  for  0.06  Hz,  the  result  shows  abou*  a  25 
degree  angle  difference  between  transformed  and  measured  DSD,  and  the  transformed 
spectral  energy  is  overestimated.  For  0.09  Hz,  the  transformed  and  measured  shallow 
waves  have  the  same  directional  spectral  distribution,  but  the  measured  data  has  greater 
energy  density  than  the  transformed.  On  the  same  day  at  0.13  Hz,  the  result  shows  good 
prediction  in  directionality.  Both  frequencies  indicate  the  transformed  DSD  have  nar¬ 
rower  distributions  compared  with  deep  water  DSD,  because  the  shinking  of  the  angle 
distribution  concentrates  the  energy  in  shallow  water  during  transformation.  The  re¬ 
fraction  effects  are  greatest  at  lower  frequency.  On  18  January  1988  for  0.09  Hz  at 
Marina,  the  results  (Figure  38)  show  a  poor  prediction,  the  directionality  having  a  95 
degree  difference,  which  is  at  low  relative  energy  within  the  spectrum  and  appears  to  be 
a  spurious  result.  On  29  January  1988,  the  results  show  a  5  degree  difference  or  less  at 
all  frequencies. 

Results  of  the  Santa  Cruz  location  are  dip-  red  in  Figure  40  -  -  The  results  ar: 
different  from  Marina.  On  18  January  1988  at  i>9  Hz,  the  predictic  as  a  15  degree 
difference,  which  is  better  than  that  at  Marina,  n  3  January  at  0.09  Hz  and  29  January' 
at  0.13  Hz  both  provide  satisfactory  results.  Table  7  summerize  the  three  days  3,  18,  and 
29  January  1988  at  different  frequencies  and  locations  comparing  transformed  and 
measured  directionality  angle  differences. 

The  variance  within  half  power  band  width  of  the  directional  spectra  at  Marina  and 
Santa  Cruz  on  all  these  days  are  given  in  Tables  8  and  9.  By  examing  Tables  7  -  9  and 
focusing  on  the  energy  peak  frequencies  for  the  three  days,  0.13  Hz  on  3  January  1988, 
0.06  Hz  on  18  January  1988  and  0.09  Hz  on  29  January  1988,  we  can  compare  the  en¬ 
ergy  and  directional  difference  between  measured  and  transformed  data  at  Marina  a:  J 
Santa  Cruz.  At  Marina,  the  directions  are  the  same  on  3  and  29  January  1988,  but  has 
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S(f.lheta),  M**2/llz/Rad 


Deep,  Shallow,  and  Transformed  DSD  at  Marina  on  1/3, 0.09  Hz 


Figure  35.  The  deep,  shallow,  and  transformed  DSD:  On  3  Jan  1988  at  Marina, 
0.09  Hz.  * 
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S(f,ilieta), 


Deep,  Shallow,  and  Transformed  DSD  at  Marins  on  1/3, 0.13  Hz 


Direction,  Deg 


Figure  36.  The  deep,  shallow,  and  transformed  DSD:  On  3  Jan  1988  at  Marina, 
0.13  Hz. 
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Iz/Rnd 


Deep,  Shallow,  and  Transformed  DSD  at  Marina  on  1/3, 0.06  Hz 


Deep,  Shallow,  and  Transformed  DSD  at  Marina  on  1/18,  0.09  Hz 


Figure  38.  The  deep,  shallow,  and  transformed  DSD:  On  18  Jan  1988  at  Marina, 
0.09  Hz.  ‘ 
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S(f.tlicta),  M**2/ll7./Rad 


Deep,  Shallow,  and  Transformed  DSD  at  Marina  on  1/29, 0.13  Hz 


Figure  39.  The  deep,  shallow, .and  transformed  DSD:  On  29  Jan  1988  at  Marina, 
0.13  Hz. 
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a  30  angle  degree  difference  on  18  January  1988,  on  the  same  day,  th  '  easur  .orgy 
is  four  times  larger  than  predicted,  but  on  3  and  29  January  1988  e  s  arc  .  same 
for  between  measured  and  predicted.  At  Santa  Cruz,  the  directional  crence  is  almost 
the  same  on  18  and  29  January  1988,  but  has  a  30  angle  degree  difference  on  3  January 
1988.  From  an  energy  point  of  view,  3  and  18  January  1988  show  the  same  energy,  but 
on  29  January  1988,  the  predicted  energy  is  almost  two  times  the  measured  energy. 

Summerizing,  the  energies  at  the  peak  frequency  were  similar  except  at  Marina  on 
the  18th  (x4)  and  Santa  Cruz  on  the  29th  (x2).  The  angles  at  the  peak  frequencies  were 
within  5  degrees  except  at  Santa  Cruz  on  the  3rd  (30°)  and  Marina  on  the  18th  (30°). 
There  is  no  obvious  correlation  with  error  on  deep  water  wave  directions  (~275°  on  the 
3rd,  300°  on  the  18th  and  245°  or.  :'ie  29th) 
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Deep,  Shallow,  and  Transformed  DSD  at  Santa  Cruz  on  1/18, 0.09  Hz 


Figure  40.  The  deep,  shallow,  and  transformed  DSD:  On  18  Jan  1988  at  Santa 
Cruz,  0.09  Hz. 


63 


S(f,t!ieia),  M**2/llz/Rad 


Deep.  Shallow,  and  Transformed  DSD  at  Santa  Cruz  on  1/3,  0.09  Hz 


Figure  41.  The  deep,  shallow,  and  transformed  DSD:  On  3  Jan  1988  at  Santa 
Cruz,  0.09  Hz. 
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V.  CONCLUSION 


Deep  water  directional  wave  spectra  were  calculated  from  data  acquired  from  the 
NDBC  buoy  and  transformed  over  the  highly  complex  bathymetry  of  Monterey  Bay  to 
shallow  water  using  a  linear  wave  refraction  model.  The  transformed  spectra  were  com¬ 
pared  with  measured  shallow  wave  spectra  at  Marina  and  Santa  Cruz. 

Two  methods  of  calculating  the  directional  spectrum  were  employed,  the  LCS  and 
EFC  techniques.  The  results  show  that  the  two  techniques  give  the  same  peak  densities 
and  energies,  but  that  the  EFC  method  has  much  superior  resolution. 

Tlv  comparison  of  measured  and  transformed  spectra  show  that  higher  frequency 
wave  c .  nditions  exhibit  more  accurate  results,  which  are  less  refracted,  and  with  rela¬ 
tively  larger  errors  occurring  at  lower  frequencies.  In  general,  linear  refraction  can  give 
reasonable  energy  and  direction  estimates  starting  with  deep  water  spectra,  but  notable 
exceptions  can  occur.  The  largest  prediction  errors  occured  at  Marina  on  18  January 
1988  for  the  case  of  a  severe  storm.  This  is  presumably  due  to  diffractive  and  non-linear 
effects  of  the  high  waves  causing  loss  of  accuracy.  The  linear  refraction  model  is  not 
suitable  for  handling  such  problems.  Further  modeling  efforts  for  Monterey  Bay  includ¬ 
ing  diffractive  and  non-linear  effects  are  needed. 
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Shallow  Water  Angle, 


APPENDIX  A.  THE  SHALLOW  WATER  ANGLE  VERSUS  DEEP 

WATER  ANGLE 


Shallow  Water  Angle  vs  Deep  Water  Angle  at  Marina,  0.09  Hz 


Deep  Water  Angle,'  Deg 


Figure  43.  The  shallow  water  angle  versus  deep  water  angle  at  Marina  0.09  IIz. 
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Shallow  Water  Angle, 


Figure  44.  The  shallow  water  angle  versus  deep  water  angle  at  Marina  0.13  IIz. 
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Shallow  Water  Angle, 


Shallow  Water  Angle  vs  Deep  Water  Angle  at  Santa  Cruz,  0.06  Hz 


Figure  45.  The  shallow  water  angle  versus  deep  water  angle  at  Santa  Cruz  0.06  IIz. 
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Shallow  Water  Angle, 


Shallow  Water  Angle  vs  Deep  Water  Angle  at  Santa  Cruz,  0.0°  Hz 
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Figure  46.  The  shallow  water  angle  versus  deep  water  angle  at  Santa  Cruz  0.09  Hz. 


Shallow  Water  Angle, 


Figure  47.  The  shallow  water  angle  versus  deep  water  angle  at  Santa  Cruz  0. 1 3  Hz. 


APPENDIX  B.  THE  DEEP,  Si  ALLOW,  AND  TRANSFORMED  DSD 


Deep,  Shallow,  and  Transformed  DSD  at  Marina  on  1/18, 0.06  Hz 


r^ure  48.  The  deep,  shallow,  and  transformed  DSD:  On  18  Jan  1988  at  Marina, 
0.06  Hz. 
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Deep,  Shallow,  and  Transformed  DSD  at  Marina  on  1/18, 0.13  Hz 


Figure  49.  The  deep,  shallow,  and  transformed  DSD:  On  18  Jan  1988  at  Marina, 
0.13  Hz. 
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S(f.theta),  M**2/Hz/Rad 
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Figure  51.  The  deep,  shallow,  and  transformed  DSD:  On  29  Jan  1988  at  Marina, 
0.09  Hz. 
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S(f,theta),  M**2/Hz/Rad 


Direction,  Deg 


Figure  53.  The  deep,  shallow,  and  transformed  DSD:  On  3  Jan  1988  at  Santa 
Cruz,  0.13  Hz. 
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Deep,  Shallow,  and  Transformed  DSD  at  Santa  Cruz  on  1/18, 0.06  Hz 


Figure  54.  The  deep,  shallow,  and  transformed  DSD:  On  18  Jan  1988  at  Santa 
Cruz,  0.06  Hz. 
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Deep,  Shallow,  and  Transformed  DSD  at  Santa  Cruz  on  1/18, 0.13  Hz 


Direction,  Deg 


Figure  55.  The  deep,  shallow,  and  transformed  DSD:  On  18  Jan  1988  at  Santa 
Cruz,  0.13  Hz. 
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S(f, theta),  M**2/Hz/Rad 


Deep,  Shallow,  and  Transformed  DSD  at  Santa  Cruz  on  1/29, 0.06  Hz 


Figure  56.  The  deep,  shallow,  and  transformed  DSD:  On  29  Jan  1988  at  Santa 
C ruz,  0.06  Hz. 
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Figure  57.  The  deep,  shallow,  and  transformed  DSD:  On  29  Jan  1988  at  Santa 
Cruz,  0.09  1  Iz. 
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